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Abstract. We discuss possible avenues to study fission dynamics starting from a 
time-dependent mean-field approach. Previous attempts to study fission dynamics 
^ I using the time-dependent Hartree-Fock (TDHF) theory are analyzed. We argue that 

(— I ' different initial conditions may be needed to describe fission dynamics depending on 

the specifics of the fission phenomenon and propose various approaches towards this 
goal. In particular, we provide preliminary calculations for studying fission following 
^ ' a heavy-ion reaction using TDHF with a density contraint. Regarding prompt muon- 

induced fission, we also suggest a new approach for combining the time-evolution of 
the muonic wave function with a microscopic treatment of fission dynamics via TDHF. 
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1. Introduction 

A fully microscopic description of the nuclear fission process is one of the most 
challenging problems in theoretical nuclear physics. Attacking this problem in 
terms of determinantal many-body wavef unctions, as is commonly done, entails the 
understanding of single-particle evolution, proper coupling to all of the relevant collective 
degrees of freedom, and incorporating the proper quantum mechanical framework. In 
principle, the semi-classical limit of this process would involve the summation over an 
ensemble of quantum mechanical fission paths that would explain the experimentally 
observed mass distributions. The problem is further exacerbated by the rich selection 
of initial conditions that lead to fission. In addition to most commonly studied forms 
of fission, namely spontaneous or neutron induced fission, which to a large degree may 
be infiuenced by the zero temperature potential energy surface, the fission events that 
occur during a heavy-ion collision, e.g. quasi-fission and fusion-fission or compound 
nucleus fission, involve excited composite systems that may not be fully equilibriated. 
For example, it has recently been shown that as we increase the temperature of the 
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system the potential energy surfaces are completely changed in comparison to the T — 
case [1]. 

The conclusion one can draw from the above discussion is that it is desirable to 
develop a dynamical approach to fission. Microscopic theories for fission have been 
discussed in the past in the context of quantized adiabatic time-dependent Hartree- 
Fock theory (ATDHF) [2] and in terms of the path-integral method [3, 4]. These 
approaches have been developed around the concept of a collective fission path and 
are biased on sub-barrier processes, mostly spontaneous fission or perhaps low-energy 
fusion. A bit more extended dynamics of fission in a collective model space has been 
studied via the time- dependent generator coordinate method employing a Hartree-Fock- 
Bogoliubov (HFB) basis with quadrupole and octupole constraints [7]. While static 
properties influencing flssion, in particular the structure of the collective potential energy 
surface, are customarily studied via adiabatic microscopic mean-fleld theories [5, 6], the 
description of excited fission channels in the classically allowed regime above the barrier 
requires truly dynamical approaches which are not confined to a collective subspace. 
The method of choice is here the TDHF theory. It becomes particularly suited for 
fission studies when coupled with a novel method of using a density-constraint to obtain 
the corresponding TDHF trajectory on the multi-dimensional potential energy surface 
of the combined nuclear system [8, 9]. Here, TDHF provides the evolution of the nuclear 
shape (density) including all of the self-consistent dynamical effects present in the mean- 
field limit. The density-constrained TDHF (DC-TDHF) approach has been successfully 
used to calculate dynamical ion-ion interaction barriers [10, 11, 12, 13, 14] for fusion 
reactions. In addition, one-body energy dissipation extracted from TDHF for low-energy 
fusion reactions was found to be in agreement with the friction coefficients based on the 
linear response theory as well as those in models where the dissipation was specifically 
adjusted to describe experiments [15]. However, while the initial conditions for a heavy- 
ion reaction are well defined, the same is not true for fission. Various types of fission 
reactions are characterized by different initial conditions which may not be obvious at 
the onset. 

In this manuscript we will first discuss previous attempts to study fission via the 
standard TDHF theory. This may provide clues to the choice of initial conditions. We 
then perform preliminary studies of fission using the DC-TDHF approach coupled with 
a collective boost. Finally, we discuss a different approach to study fission dynamics, 
namely prompt muon-induced fission involving nuclear dynamics described by TDHF. 

2. Previous TDHF Studies and Initial Conditions 

In general there are very few attempts to generalize the static mean-field approach to 
a dynamical approach for fission. In this section we will discuss some of the previous 
attempts to study fission using the standard TDHF approach. 

Perhaps the most well known dynamical fission study is the one discussed in 
Ref. [16]. Here, the authors studied the fission of ^^^U using TDHF with axial and 
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reflection symmetries and with a simplified form of the Skyrme interaction without 
the spin-orbit force. The initial state was determined via a constrained Hartree-Fock 
(CHF) calculation to be 1 MeV beyond the saddle point. Furthermore, due to the 
imposed axial symmetry restriction a time-dependent pairing had to be introduced in 
order to force the coupling of otherwise unmixed states. This is an important point 
to stress because the significance of pairing for the general fission problem is still an 
open question. Here, pairing was solely introduced to break unphysical symmetries. As 
anticipated, the authors do find a number of fission paths that strongly depend on the 
choice of the pairing properties. This seems to suggest that to study dynamics of fission 
no symmetry restrictions should be made. However, this study does not directly address 
the problem of how the system evolves to the chosen initial state and naturally suffers 
from too many approximations made to simplify the numerical calculations. 

Another early study of fission via TDHF is given in Ref. [17]. Here, the authors 
perform a similar study except the initial CHF states are chosen from single-center or 
two-center determinantal trial functions and no pairing. These states are then boosted 
with a collective quadrupole velocity field. One observes that the two-center initial 
states can be induced to fission while the single-center calculations do not result in 
fission despite the deposition of large amount of collective energy. The authors draw the 
conclusion that the problem lies with having a single Slater determinant and the lack 
of understanding as to which degree of freedom to invoke for the correct initialization 
of fission dynamics. 

Finally, an earlier TDHF study limited to slabs of nuclear matter [18] and using 
a velocity field to trigger fission leads to the observation that starting from the HF 
ground-state does not produce fission for reasonable velocity fields, while starting from 
HF states that are obtained by hfting some nucleons from occupied states to previously 
unoccupied levels fission can be invoked with relative ease. 

The conclusion one can draw from all these studies is that the initial many-body 
state for fission should be some form of a correlated state that contains couplings to 
continuum or other excited states. The actual construction of this initial state will 
also depend on the type of fission that one is planning to study, prompt fission being 
perhaps the most difficult one since here the correlations and virtual excitations are 
the sole drivers to fission. For the case of neutron induced fission, the selection of 
where the neutron energy is deposited in terms of the single-particle picture and the 
construction of the resulting correlated state is required. The case of fusion-fission and 
especially quasi-fission can be better understood in terms of the states formed during 
a TDHF collision since these do contain multiple centers and large couplings to the 
continuum. This was noted by a statement attributed to A. K. Kerman in Ref. [19] 
as: ...it is tempting to speculate that states formed in TDHF fusion reactions might be 
associated with single-particle wave functions which are approximate eigenstates of the 
instantaneous HF Hamiltonian and the many-body wave function constructed from these 
single-particle wave functions could then be considered a transition state to processes that 
are not taken into account in TDHF theory. The problem these authors faced was the 
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unavailability of a numerical method to extract the undissipated energy from TDHF 
dctcrminantal wave functions since evolving to infinite time was not an option. They 
tried the imaginary-time method to find the corresponding static solutions but the 
numerical procedure always lead to the ground state of the combined system. This now 
can be accomplished using the density constraint combined with the imaginary-time 
method, as it used in DC-TDHF calculations and will be discussed in the next section. 

3. TDHF with Density Constraint 

In order to validate some of the above findings we have tried to achieve dynamical fission 
by starting from either the ground state or the traditional fission isomer obtained by 
CHF calculations with a quadrupole constraint. Specifically, we have looked at ^^"^Pu 
and ^^^U systems. The difference from earlier calculations was that here fully three- 
dimensional TDHF codes with no symmetry assumptions and the full modern Skyrme 
forces were used. Using collective boost operators of the type 6^'^^'^^°^^^ we were unable to 
obtain fission for reasonable excitation energies. If a large amount of collective energy is 
deposited we observe symmetric fission, which for the case of ^^°Pu is the least probable 
outcome. Wc arc still investigating this approach by using a combination of quadrupole 
and octupolc boost operators in an attempt to channel the system to the appropriate 
fission path. However, it is becoming apparent that the energy may have to be deposited 
in a more selective way than to the entire nucleus. Wc arc also studying the possibility 
of constructing excited many-body states and using them as an initial state for fission. 
Naturally, all of these studies only apply to the case of neutron induced fission. 

As we have also discussed above, fission following a heavy-ion collision (quasi- 
fission or fusion-fission) may be studied by creating the initial states via the actual 
TDHF evolution of the collision followed by density constraint calculation. Wc have 
first performed collisions of ^^'^Zr+^^'^Xe, which results in the extensively studied ^"^^Pu 
composite system. The choice of these nuclei were motivated by the most probable 
fragment mass and charge for the induced fission of ^^''Pu. Along the TDHF trajectory 
the density-constraint was applied to obtain the ion-ion potentials. In Fig. la we plot the 
result at three different bombarding energies. As wc sec at the lowest energy the system 
undergoes a deep-inelastic collision, while at the highest energy we observe complete 
fusion, as seen from the fact that the system is void of any collective motion for a long 
time (2500 fm/c). On the other hand, for the intermediate energy the system executes 
large collective oscillations with very little damping for the same time interval. This 
is what we call the transition state and indicates that the system is trapped in an 
isomeric minimum and executes collective motion. Following this idea we have started 
from a density along this path and performed an unconstrained minimization. The 
result was an isomer of the ^''^Pu system shown in Fig lb (/32 = 2.27, Q20 = 2306, and 
Q30 ~ —286^/^). This isomer is far from the well known isomer of ^^°Pu, which cannot 
be reached in such reactions. One feature of these transition states is that they fission 
with a very small collective boost, e.g. e^'^P'^^°^^\ Using p — 0.0025 the nucleus fissions 
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Figure 1. 

(a) Ion-ion potential for the ^'^"Zr+^'*°Xe system at three different energies. 

(b) New shape isomer of ^*°Pu fomid from the transition state. 



with an initial excitation energy of 7 MeV. By using the above boost and performing 
DC-TDHF calculations for the fissioning nuclei one can investigate the potential barriers 
in the vicinity of the isomeric minimum. This is shown in Fig. 2. The resulting masses of 
the two fragments are A\^Z\ = 106, 42 and A2, Z2 = 134, 52. We feel that this approach 
to study fission following a heavy-ion reaction is very promising and further explorations 
will be done in the near future. 
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Figure 2. Potential barrier in the vicinity of the new shape isomer. 



4. Fission "clocks" 

In this section we will discuss a different approach to study fission dynamics, namely 
muon-induced fission, which may also shed light on to nuclear energy dissipation defined 
as the irreversible flow of energy from collective to intrinsic single-particle motion [20]. 
There are two different mechanisms that contribute to nuclear energy dissipation: two- 
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body collisions and one-body friction. The latter mechanism is caused by the moving 
walls of the one-body time-dependent mean field as described by TDHF. The role 
played by these two dissipation mechanisms in fission and in heavy-ion reactions is 
not yet completely understood. Assuming that friction is caused by two-body collisions, 
Davies et al. [21] calculated the effect of viscosity on the dynamics of fission; they 
extracted a viscosity coefficient fj, — 0.015 Tera Poise from a comparison of theoretical 
and experimental values for the kinetic energies of fission fragments. The corresponding 
time delay for the nuclear motion from the saddle to the scission point was found to be 
of the order of At = 1 x 10~^^ s. However, in one-body dissipation models the time 
delay is an order of magnitude larger. 

Several experimental techniques are sensitive to the energy dissipation in nuclear 
fission. At high excitation energy, the multiplicity of pre-scission neutrons [22] or 
photons [23] depends on the dissipation strength. At low excitation energy, the process 
of prompt muon-induced fission [24, 25, 26, 27, 28, 29] provides a suitable clock. 

Following the formation of an excited muonic atom, inner shell transitions may 
proceed by inverse internal conversion where the muonic excitation energy is transferred 
to the nucleus. In actinides, the 2p — > Is and the 3d — >■ Is muonic transitions result 
in excitation of the nuclear giant dipole and giant quadrupole resonances, respectively, 
which act as doorway states for fission. The nuclear excitation energy is typically 6.5 — 10 
MeV, i.e. E* exceeds the fission barrier. Because the muon lifetime is long compared 
to the timescale of prompt nuclear fission, the motion of the muon in the Coulomb 
field of the fissioning nucleus may be utilized to learn about the dynamics of fission. 
If there is large friction between the outer fission barrier and the scission point the 
muon will remain in the lowest molecular energy level and emerge in the Is bound state 
of the heavy fission fragment. On the other hand, if friction is small, there is a non- 
zero probability that the muon may be promoted to higher-lying molecular orbitals from 
where it will end up attached to the hght fission fragment. Therefore, theoretical studies 
of the muon-attachment probability to the light fission fragment, Pl, in combination 
with experimental data can be utilized to analyze the dynamics of fission. 

From a theoretical point of view, prompt muon-induced fission has several attractive 
features. Because E* exceeds the fission barrier we do not face the difficult problem of 
a microscopic description of barrier tunneling. The muon dynamics is determined by 
solving either the Schrodinger equation [25, 26] or the Dirac equation [28, 29] in the 
presence of a time-dependent external Coulomb field which is generated by the nuclear 
charge density during fission. So far, the motion of the fissioning nucleus has been 
described by a simple classical model involving a linear friction force. 

Our time- dependent Dirac equation calculations [29] predict a strong mass 
asymmetry dependence of the muon attachment probability to the light fission 
fragment, in agreement with experimental data. The theory also predicts a dependence 
of Pl on the dissipated energy. By comparing our theoretical results to the experimental 
data of ref. [27] we extract a dissipated energy of order — 10 MeV for ^^'''Np. The 10 
MeV value corresponds to a fission time delay from saddle to scission of order 2 x 10~^^ 
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s. This value of E^iss = 10 MeV agrees with results from other low-energy fission 
measurements that are based on the odd-even effect in the charge yields of fission 
fragments [30]. 

In the future, we might be able to test the validity of one-body (mean-field) energy 
dissipation in large- amplitude collective motion by combining the time-evolution of the 
muonic wave function (as described by the 3-D Dirac equation) with a microscopic 
treatment of fission dynamics with our 3-D TDHF code. 

5. Summary 

The study of fission dynamics represents one of the most challenging problems in nuclear 
many-body theory. In the absence of a full microscopic and quantum mechanical theory 
for the many-body tunneling process we ask the question about the generalization 
of the mean-field based calculations, that have been extensively used to calculate 
potential energy surfaces and fission barriers, to the dynamical case using a TDHF based 
approach. A careful analysis leads to the conclusion that this may be possible provided 
the proper initial conditions for the dynamical calculations are chosen. We show that 
the most straightforward study using TDHF may be fission which follows a heavy-ion 
reaction. This could be achieved via the use of density constraint and collective boost 
operators on top of the TDHF dynamics. Studying the dynamics of fission induced by 
low energy neutrons requires the construction of the initial state which represents the 
state of the nucleus after the energy transfer. One could construct such excited states 
by using standard variational procedures, which we are currently investigating. Clearly, 
all of these initial states will contain correlations beyond the simple CHF state. Exactly 
what type of correlations are needed could emerge from such studies. 

Acknowledgments 

This work has been supported by the U.S. Department of Energy under grant No. DE- 
FG02-96ER40963 with Vanderbilt University, and by the German BMBF under contract 
Nos. 06FY159D and 06ER142D. 

References 

[1] Pei J C, Nazarewicz W, Sheikh J A and Kerman A K 2009 Phys. Rev. Lett. 102 192501 
[2] Goeke K, Griimmer F and Reinhard P-G 1983 Ann. Phys. 150 504 

[3] Negele J W 1982 Rev. Mod. Phys. 54 913 

[4] Gocke K, Reinhard P-G and Reinhardt H 1982 Nucl. Phys. A 378 474 
[5] Dubray N, Goutte H and Delaroche J-P 2008 Phys. Rev. C 77 014310 
[6] Schindzielorz N et al. 2009 Int. J. Mod. Phys. E 18 773 

[7] Goutte H, Berger J F, Casoh P and Gogny D 2008 Phys. Rev. C 71 024316 

[8] Cusson R Y, Reinhard P-G, Strayer M R, Maruhn J A and Greiner W 1985 Z. Phys. A 320, 475 
[9] Umar A S, Strayer M R, Cusson R Y, Reinhard P-G and Bromley D A Phys. Rev. C 32 172 
[10] Umar A S and Oberacker V E 2006 Phys. Rev. C 74 021601 (R) 



Microscopic Fission Dynamics 



8 



[11] Umar A S and Obcracker V E 2006 Phys. Rev. C 74 061601(R) 

[12] Umar A S and Oberacker V E 2007 Phys. Rev. C 76 014614 

[13] Umar A S and Oberacker V E 2008 Phys. Rev. C 77 064605 

[14] Umar A S and Oberacker V E 2009 Eur. Phys. J. A 39 243 

[15] Washiyama Kouhci, Lacroix Denis and Ayik Sakir 2009 Phys. Rev. C 79 024609 

[16] Negele J W, Koonin S E, MoUer P, Nix J R and Sierk A J 1978 Phys. Rev. C 17 1098 

[17] Okolowicz J, Irvine J M and Nemeth J 1983 J. Phys. G 9 1385 

[18] Dietrich K and Ncmcth J 1981 Z. Phys. A 300 183 

[19] Davies K T R, Sandhya Devi K R and Strayer M R 1981 Phys. Rev. C 24 2576 
[20] Hasse R W 1978 Rep. Prog. Phys. 41 1027 

[21] Davies K T R, Sierk A J and Nix J R 1976 Phys. Rev. C 13 2385 

[22] Gavron A et al. 1987 Phys. Rev. C 35, 579 
[23] Hofmann D J ei al. 1995 Phys. Rev. C 51 2597 

[24] Wilcke W W, Johnson M W, Schroder W U, Huizenga J R and Perry D G 1978 Phys. Rev. C 18 

1452 

[25] Maruhn J A, Oberacker V E and Maruhn-Rezwani V 1980 Phys. Rev. Lett. 44 1576 
[26] Glanders P, Nilsson S G and Moller P 1980 Phys. Lett. B 90 193 
[27] Risse F et al. 1991 Z. Phys. A339 427 

[28] Oberacker V E, Umar A S, Wells J C, Bottcher C, Strayer M R and Maruhn J A 1993 Phys. Rev. 
C 48 1297 

[29] Oberacker V E 1999 Heavy-Ion Physics Vol. 10, No.2-3, 221 

[30] Wagemans C 1991 The Nuclear Fission Process (CRC Press, Boca Raton) p. 418 



